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Abstract 

We analyze the resonant inelastic x-ray scattering (RIXS) spectra at the Ni K edge in an anti- 
ferromagnetic insulator NiO by applying the theory developed by the present authors. It is based 
on the Keldysh Green's function formalism, and treats the core-hole potential in the intermediate 
state within the Born approximation. We calculate the single-particle energy bands within the 
Hartree-Fock approximation on the basis of the multi-orbital tight-binding model. Using these 
energy bands together with the 4p density of states from an ab initio band structure calculation, 
we calculate the RIXS intensities as a function of energy loss. By taking account of electron corre- 
lation within the random phase approximation (RPA), we obtain quantitative agreement with the 
experimental RIXS spectra, which consist of prominent two peaks around 5 eV and 8 eV, and the 
former shows considerable dispersion while the latter shows no dispersion. We interpret the peaks 
as a result of a band-to-band transition augmented by the RPA correlation. 

PACS numbers: 78.70.Ck 71.20.Be 71.28.+d 78.20.Bh 
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I. INTRODUCTION 



Excitations in solids are fundamental to describe physical properties such as the response 
to external perturbations and temperature dependence. They may be characterized into 
two types, spin and charge excitations. For the former, the inelastic neutron scattering is 
quite powerful to investigate energy-momentum relations. By contrast, charge excitations 
have been investigated by measuring the optical conductivity, but the momentum transfer 
is limited to nearly zero.- The electron energy loss spectroscopy can detect the momen- 
tum dependence of charge excitations, but it suffers from strong multiple scattering effect s.- 
Recently, taking advantage of strong synchrotron sources, the resonant inelastic x-ray scat- 
tering (RIXS) has become a powerful tool to probe charge excitations in solids.-^^ 1 ^ 
In transition-metal compounds, i^-edge resonances are widely used to observe momentum 
dependence, because corresponding x-rays have wavelengths the same order of lattice spac- 
ing. The process is described as a second-order optical process that the Is core electron is 
prompted to an empty 4p state by absorbing photon, then charge excitations are created in 
order to screen the core-hole potential, and finally the photo-excited 4p electron is recom- 
bined with the core hole by emitting photon. In the end, charge excitations are left with 
energy and momentum transferred from photon. In cuprates, the RIXS spectra are found 
to have clear momentum dependence.-^ 

For NiO, an pioneering Ni K-edge RIXS experiment has been carried out by Kao et al., 
in which the spectral peaks were not well resolved had no clear momentum dependence, 
probably due to the experimental resolution.- Recently, a new K-edge RIXS experiment 
has been carried out, at Taiwan beam line in SPring-8.— They have observed the spectra 
as a function of energy loss by tuning the incident photon energy at 8351 eV, which corre- 
sponds to the Ni i^-edge absorption peak. The observed spectra consist of two prominent 
peaks at 5 eV and 8 eV, and the 5-eV-peak considerably changes while the 8-eV-peak does 
not change with changing momenta. In addition, extra tiny peaks are found below 4 eV, 
which are called u d-d excitation" . In this paper, we analyze the RIXS spectra for NiO by 
developing the formalism of Nomura and Igarashi.— ^ This theory is based on the many- 
body formalism of Keldysh, and is regarded as an extension of the resonant Raman theory 
developed by Nozieres and Abrahams.— The core-hole potential is treated within the Born 
approximation. Higher-order effects beyond the Born approximation have been evaluated 
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on the K-edge RIXS in L^CuO^— Although the core-hole potential is rather strong, the 
higher-order effects are found to cause only minor change in the spectral shape.— In this 
situation, the RIXS spectra can be connected to the 3d-density-density correlation function 
in the equilibrium state. We develop the formalism by clarifying the equivalence of the 
Keldysh formalism and the conventional Green's function formalism and by deriving the 
RIXS formula for possible bound-states corresponding to the d-d excitations. Advantages 
of the present formalism are, in contrast to the numerical diagonalization method, that it 
is applicable to three-dimensional models consisting of many orbitals, and that it provides 
clear physical interpretation. 

We construct a detailed multi-orbital tight-binding model including all 3d orbitals as well 
as the full intra-atomic Coulomb interaction between 3c? orbitals. Applying the Hartree-Fock 
approximation (HFA) to the model, we obtain the antiferromagnetic (AF) solution with an 
appropriate description of the single-particle spectra having a large energy gap. Unoccupied 
3d states on a Ni site consist of minority spin states on the site and have almost the e 9 - 
character. Note that the band calculation with the local density approximation (LDA) 
fails to reproduce the energy gap of 4 eV. The RIXS spectra are interpreted as a result 
of a band-to-band transition to screen the Is core hole. Therefore the transition occurs 
through the amplitude of the e 9 -character. Treating electron correlations by the random 
phase approximation (RPA), we find that spectral shape as a function of energy loss is 
strongly modified in the continuum spectra. In order to obtain quantitative agreement with 
the experiment, we need to take account of the RPA correlation. Note that the same analysis 
of the RIXS in cuprates has been successful to reproduce the experimental spectra. The 
success of the analysis for NiO would add another evidence of the usefulness of the present 
scheme of analyzing the RIXS spectra. 

The present formalism describes the "d-d" excitation by the bound state in the density- 
density correlation function. We obtain bound states but with extremely small intensities. 
We will discuss possible reasons for the discrepancy with the experiment. 

The present paper is organized as follows. In Sec. UH we study the electronic structure 
on the basis of the band calculation, and then introduce a multi-orbital tight-binding model. 
We calculate the single-particle Green's function within the HFA in the AFM phase of NiO. 
In Sec. [Till we summarize the formula for the RIXS spectra, including the discussion of 
the bound states. In Sec. ITVl we calculate the RIXS spectra by taking account of the 
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Figure 1: Schematic view of a NiO crystal with type-II AF order. Only Ni atoms are shown. Wave 
vector Qaf characterizes the AF modulation. The direction of magnetic moment on the Ni site 
denoted by filled circle is antiparallel to that on the site denoted by open circle. 

RPA correction in comparison with the experiment. Section |V] is devoted to the concluding 
remarks. 

II. ELECTRONIC STRUCTURE OF NICKEL OXIDE 

The crystal structure of NiO is the NaCl-type with the lattice constant of a = 4.177 
A. Ni atoms form an fee lattice, as shown in Fig. [TJ Type -II AF order develops below 
Tn = 523 K. The order parameter is characterized by a wave vector directing to one of four 
body diagonals in the fee lattice.— 

A. Ab initio calculation 

We calculate the electronic band structure using the muffin-tin KKR method within 
the LDA. Although we obtain a stable AF self-consistent solution, the energy gap is much 
smaller than the experimental one, ~ 4.3 eV.— Figure [2] shows the calculated density of 
states (DOS) projected onto Ni3<i and 02p states. The difference between the calculation 
and the experiment in the energy gap may be improved by using the so called LDA+[7 
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Figure 2: Density of states projected onto Ni3d and 02p states calculated within the LDA. The 
energy zero is at the top of the valence band. 

method, which is a hybrid of the LDA and the HFA for the Coulomb interaction in the Ni 
3d orbitals. Since we need to calculate the two-particle correlation function within the RPA, 
we introduce a multi-orbital tight-binding model in place of the ab initio calculation, and 
apply the RPA in the following. 

In contrast to the failure for the 3d states, we expect that the 4p bands are well described 
within the LDA, since the 4p band has a wide width > 20 eV and thereby are weakly 
correlated. We show the 4pDOS convoluted with the Lorentzian function with the full 
width of half maximum (FWHM) 2 eV in Fig. [3j The FWHM corresponds to the core-hole 
life-time width. The calculated curve agrees fairly well with the experimental one. 

B. Multi-orbital tight-binding model 

We introduce a multi-orbital tight-binding model defined by 



H = H + Hi, 



+ E E {fftJvV^Vj'l'a + tf-C.) + E E (^i'm'dlmA'm'a + H.C 
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Figure 3: Density of states projected onto Ni4p states calculated within the LDA. The experimental 
curve is taken from Fig. 4 in ref. 1 171 . 

The part H represents the kinetic energy, where d ima and p^ a denote the annihilation 
operators of an electron with spin a in the 3d orbit m of Ni site i and the annihilation operator 
of an electron with spin a in the 2p orbit I of the O-site j, respectively. Number operators 
nf ma and n p jia are defined by nf ma = d\ ma d ima , n p jl<T = p} ea p jerT - The transfer integrals, t*^, 
tj£ji£i tim,i'm'i are evaluated from the Slater-Koster two-center integrals, (pda), (pdir), (ppcr), 
(pprc), (ddcr), (ddn), (dd5). 18 The (i-level position relative to the p-levels is given by the 
charge-transfer energy A defined by A = — E p + 8U for the d 8 configuration.— Here U is 
the multiplet-averaged d-d Coulomb interaction given by U = F° — (2/63) F 2 — (2/63) F 4 , 
where F°, F 2 , and F 4 are Slater integrals for 3d orbitals. The part Hi represents the intra- 
atomic Coulomb interaction on TM sites. The Coulomb interaction on O sites is neglected. 
The interaction matrix element g iv\Vi] ^3^4) is written in terms of F°, F 2 , and F 4 [y stands 
for (m, a)). 

We determine most parameter values from a cluster-model analysis of photo-emission 
spectra.— The values for (dda), (ddir), and (dd5) can not be determined from the cluster- 
model analysis, and therefore we set them close to Mattheiss' LDA estimates.— Among 
Slater integrals, F 2 and F 4 are known to be slightly screened by solid-state effects, so that 
we use the values multiplying 0.8 to atomic values. On the other hand, F° is known to 
be considerably screened, so that we regard the value as an adjustable parameter to get a 
reasonable band gap. Table [I] lists the parameter values used in the present calculation. 
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Table I: Parameter values for the tight-binding model of NiO in units of eV. 





SK 


par am. 


Slater Integral 


Ni 


dda 


-0.227 


F° 5.00 




ddir 


-0.103 


F 2 10.57 




dd5 


-0.010 


F 4 7.56 


Ni-0 


pdir 


-1.400 


charge transfer energy 




pda 


0.630 


A 4.5 





ppa 


0.600 






pp-K 


-0.150 





C. Hartree-Fock Approximation 

For the AFM order shown in Fig. [T], a unit cell contains two Ni atoms and two O atoms. 
Labeling a unit cell by i], we introduce the Fourier transform of the annihilation operator 
Arm ( = or Vjla) i n the magnetic Brillouin zone, 

V 

Here represents a position vector of the unit cell, and rj runs over N/2 unit cells. Note that 
a single phase factor kr^ is applied to all the states in each unit cell. With these operators, 
the single-particle Green's function is introduced in a matrix form, 

[6(k,w)U = -if (T(A k ,(t)A^,(0)))e^dt (6) 

In the HFA, we disregard the fluctuation terms in Hj and approximate Hj by 

H ? F = \Y. E r(°)( i , 1 ^^3^)(4 2 ^a)4 1 ^4, (7) 

where T^is the antisymmetric vertex function defined by 

r (0) (z/!Z/ 2 ; z/ 3 ^ 4 ) = g(v\V2\ Wi) - g{viV2, u 4 u 3 ). (8) 

where (A) denotes the ground state average of the operator A. With the Hamiltonian 
+ Hf F , we obtain the following relation: 

(ujI - J(k))G(k, u) = I, (9) 
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Figure 4: Energy band as a function of momentum along symmetry lines within the HFA. The 
energy zero is at the top of the valence band. 

where I is the unit matrix, and J(k) is given by 

[H + Hf F , A k/ J = -V( k Mv- (10) 

Introducing an unitary matrix £/(k) to diagonalize J(k), that is, [?7(k) _1 J(\s)U(\c)]jji = 
Ej(k)5jj', we express the Green's function as 

G(k, u) = U{k)D(k, u)U{k)-\ (11) 

with 

[£){K u)]i >' = u-EM + rngziEiM-to) 6 "'' (12) 

where /xo is the chemical potential. The J(k) contains expectation values on the ground 
state, which should be self-consistently determined from the relation, 

«<V> = f E / [G(k,o;)]^ 0+ ^, (13) 

where £ and £' stand for /i and /i', respectively, at one of two Ni sites in a unit cell. Only 
the diagonal parts are non-vanishing, if 3d orbitals are specified by m = xy, yz, zx, x 2 — y 2 , 
3z 2 — r 2 with x, y, z referring to the crystal axes, that is, (d iv di V i) ^ only for v = v' . 

A stable self-consistent solution exists for the AFM order shown in Fig. [H Figure H] 
shows the energy band as a function of momentum along symmetry lines. The energy gap 
is obtained as ~ 4 eV, in consistent with the experiment. 
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Figure 5: Spin-resolved DOS projected onto Ni 3d states and onto O 2p states. The energy zero is 
at the top of the valence band. 

Figures and E] shows the DOS projected onto Ni 3d and O 2p states, and the 3dDOS 
divided into the e g and t 2g characters, respectively. The states around the top of the valence 
band have relatively large weights of O 2p state, implying that NiO is an insulator of charge 
transfer type. Regarding the majority spin states, the 3d states are almost fully occupied; 
the e g character is concentrated around the top and bottom of the valence band, while the 
t 2g character is around the middle of the valence band. Regarding the minority spin states, 
the unoccupied states have almost the e g character. The t 2g character dominates around the 
top of the valence band, while the e g character are widely distributed with relative small 
weight. As become clear later, the electron and the hole created in the RIXS process have 
the same spin. Therefore, the pair creation in the RIXS takes place in the minority spin 
states with the e g character. 

The electron correlation modifies the bands given by the HFA. One of the most prominent 
difference is that a "satellite" peak is created around 9 eV below the top of the valence 
band. In addition, the 3d states are pushed to upper energy position as a counter effect of 
the satellite creation. However, these modifications are limited deep in the low energy part 
of the valence band, mostly with the majority spin. Therefore, the electron correlation on 
the single-particle excitation may have little influence on the RIXS spectra. 
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Figure 6: Spin-resolved DOS divided into the t2 g and e g characters. The energy zero is at the top 
of the valence band. 

III. FORMULA FOR RIXS SPECTRA 
A. general expression 

For the interaction between photon and matter, we consider the dipole transition at the 
K edge, where the ls-core electron is excited to the 4p band with absorbing photon and the 
reverse process takes place. This process may be described by 



where &q represents the 77-th component (77 = x, y, z) of two kinds of polarization vectors 
(a = 1,2) of photon. Since the Is state is so localized that the Is — > Ap dipole transition 
matrix element is well approximated as a constant w. Annihilation operators p'^ and Sj a 
are for states Ap v and state Is at Cu site j, respectively. The Hamiltonians for the core 
electron and for 4p electrons are given by 




(14) 




(15) 




(16) 



kr/cr 
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The photo-created ls-core hole induces charge excitations through the attractive core-hole 
potential, which may be described by 

H ls - 3 d = V 4m* d ima4<T' S i<r'- ( 17 ) 

imaa' 

Here V may be 5 — 10 eV in NiO. 

This process is diagrammatically represented in Fig. [71 where the Born approximation is 
applied to the screening of the core-hole potential. The shaded part in the figure represents 
the Keldysh-type Green's function, 

Y+s(d,s'-s) = I r+ c -(q,u,)e-M*'-)^ (is) 

= (PqmV'(s>qL(s)), (19) 

where superscripts + and — stand for the backward and outward time legs, respectively- 
Confining our discussion to the zero temperature, we take the average over the ground state. 
Therefore, Y + ~ is nothing but a conventional correlation function in the equilibrium state. 
The superfix A on the density operator represents A or B sites of Ni in a unit cell; we have 
for \ = A, 

P%na — j ^ ^k+qmcr ^kmcr • (20) 
k 

The momentum conservation requires the relation q = q» — q/, and k runs over the magnetic 
first BZ. 

The product of Green's functions for the 4p electron and for the core hole, which are rep- 
resented by lines with labels "4p" and "Is" in the figure, gives simply a factor exp[i(el p (p) — 
€xs — %£\ s — Wi)t] on the outward time leg and a factor exp[— i(eL(p) — eia+«Ti a — uJi)(t' — u')) 
on the backward time leg, where T ls is a lifetime broadening width of the Is core hole. The 
Keldysh Green's function carries the time dependent factor e luJS to the outward time leg 
and e~ luls ' to the backward time leg. Note that the core-hole potential works only in inter- 
vals [t,0] and Integrating the time factor combined to the above product of Green's 
functions, with respect to s and t in the region of t < s < 0, — oo < t < 0, we obtain 

/0 pO 
dt ex P[*(4,(P) - ei- - - "i)t] / e l " s 
■oo „ Jt 

Vpl p (e)de 



p 



(Ui + e ls + iTi s - e)(ui - u + e ls + iT u - e) ' 



(21) 
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Figure 7: Diagrams for the RIXS intensity within the Born approximation to the Is core-hole 
potential. The dotted lines represent the core-hole potential V . The shaded part represents the 
Keldysh-type Green's function, which connects the outward time leg on the top half and the 
backward time leg on the bottom half. 

Here the sum over 4p states is replaced by the integration of the Ap DOS projected onto the 
7] (= x, y, z) symmetry, pj p (e). A similar factor has been derived in third-order perturbation 
theory by Abbamonte et al.— The integration with respect to s' and t' in the backward 
time leg gives the term complex-conjugate to Eq. (T2TT) . The integration with respect to 
u' give the energy conservation factor, which guarantees that to in Eq. (12T1) is the energy 
loss, uj = io i — ujf. Finally, combining these relations together, we obtain an expression of 
the RIXS intensity for the incident photon with the momentum and energy q% = (q^cuj), 
polarization e^, and the scattered photon with the momentum and energy qj = (q^,a;/), 
polarization e^: 

2 

(22) 

Here q = (q, uj) with q = q^ — q/, oj = Ui — LOf. Suffix £ and £' represent the 3d states 
distinguishing the A or B Ni sites. 
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W{q h an\ q f , at f ) = V Y+ (q) 
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B. The correlation function 



We consider diagrams where one electron-hole pair remains after the RIXS process, as 
shown in Fig. E](a). We write Y^j(q) in the form, 

Y &® = E *i*a«)K*IuW^m, (23) 

6666 

with 

= 27r^^^tf(w-E i ;(k + q)+^ i (k))[l-^(k + q)]fi i (k) 

x £^(k + q)E&,(k + q)^ 4 ,(k)^(k), (24) 

where j and j' denote energy eigenstates. The ^-function in n + ~(°)(g) indicates that the 
inter band transition from the valence band to the conduction band gives rise to the RIXS 
intensity. Only the weight of 3c? states in the bands contributes to the intensity. We evaluate 
the vertex function in Eq. f[2"3"j) within the RPA. Collecting the ladder diagrams shown in 
Fig. Mjo), we obtain 

A66,S(?)= \i-tF--(q)]~ 1 (25) 
where F (q) is the conventional time-ordered propagator given by 

[^~(<?)ki6,66 = / ^=- G uU^ k o) G ^( k + ^ k o + UJ ) 

k J 

2 J2 ^(k)^(k)%^(k + q)L&,(k + q) 



N 

k 



rij (k) [1 - rif (k + q)] rij, (k + q) [1 - n, (k)] 



.(26) 



u- E j/ (k + q)+E j (k)+iS lu - E r {k + q) + ^(k) - i5 _ 

The four-point vertex is non-zero only for £ 1; £ 2 > £ 3 , £ 4 belonging to the same Ni site, that 
is. 

[r]66,66 = r (0) (£i£4;66), (27) 

and otherwise it is zero. 

As already pointed out, the Keldysh-type Green's function Y^T(g) is equivalent to the 
conventional correlation function. Therefore, it may be more convenient to analyze the 
function with the help of the conventional time-ordered Green's function, 

*&(q,t) = -i(T[p* m ,AtMm<r)HO)}), (28) 
13 
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Figure 8: (a)Diagrams for (q). Solid lines with arrows are not the time-ordered Green's 
functions but the Keldysh-type ones connecting the outward time leg with the backward time leg. 

(b) Vertex function A(q) within the RPA. Solid lines with arrows are the conventional time-ordered 
Green's functions. The square represents the four-point vertex of the 3d Coulomb interaction 
effective only at the same Ni sites. (c)Diagrams for the time-ordered Y T (q) within the RPA. 

with T being the time-ordering operator. Considering the diagrams shown in Figs. 0(b) and 

(c) , it is expressed within the RPA as 

= {f-(q)[i - tp ~~W%e# ■ ^ 

This expression leads to Eq. fl2"31) with the help of the fluctuation-dissipation theorem (FDT), 

To show this fact, we first rewrite Eq. (|2U|) as 

F-(q)=Fr{q)+iFr(q), (31) 

where F{ (q) and _F 2 (g) are Hermitian matrices. Then, using the Hermitian property, we 
transform Y^(q)* as 

Y^(qT = {[I-(Fr-(q)-iF 2 -(q))r]- 1 [F l -(q)-iF 2 -(q)] 

x [I- T(Fr(q) + iF 2 -(q)W ~ WO?) + i^TC?))]" 1 }^- (32) 

Combining the similar expression for Y^,(g), we have 

(33) 
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Since — 2F 2 (g) is equivalent to U^°\q) for uj > 0, the RXS intensity given by Eq. ( 1331) 
through Eq. (1301) is equivalent to the one given by Eq. ( 1231) . 

Equation (1291) has sometimes poles for some frequencies below the energy continuum of 
an electron-hole pair. These bound states may be called as u d-d excitations" , and give rise 
to extra RIXS peaks. To analyze the problem, we rewrite Eq. (1291) as 

(?) ' r . i ,', ii - (34) 

where F (q)^ 1 becomes an Hermitian matrix for u below the energy continuum. Therefore, 
F (q)^ 1 — r can be diagonalized by a unitary matrix. Let the diagonalized matrix have 
a zero component at u = Wb((1) with the corresponding eigenvector B^ 2 (q). Then, Y^(q) 
can be expanded around to ~ ^(q) as 

*S<9> = -^4- v (35) 

with 

B^r(q)BUq) 

Q's(q) = ! • (36) 



Substituting the right hand side of Eq. fl30l) by this equation, we obtain 

^2 Y+f (g) = 2tt ) % - w fl (q) ) . (37) 

This relation is to be inserted into Eq. (1221 for evaluating the RIXS spectra. 



IV. CALCULATED RESULTS 

We calculate the RIXS intensity from Eq. (122]) . The 5-function in Y + ~(q) is replaced 
by the Lorentzian function with the FWHM = 0.2 eV in order to take account of the 
instrumental resolution. Corresponding to the experimental situation that the polarization 
of photon is not specified, we use the realistic 4p-DOS in evaluating L v B (uJi] u) (Eq. ( J2~T1) ) and 
tune the incident photon energy on the top of the 4p-DOS. Figure [9] shows the calculated 
spectra as a function of energy loss in comparison with the experiment by Ishii et al.~° The 
momentum transfer is chosen along symmetry axes in the Brillouin zone. 
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Energy Loss (eV) Energy Loss (eV) Energy Loss (eV) 

Figure 9: RIXS spectra as a function of energy loss u = Ui — tUf, in comparison with the experiment 
by Ishii et al.— Momentum transfer is along the symmetry lines: (a) T- ^(0.25, 0.25, 0.25)-L, (b) 
r-^(0.3,0,0)-X, and (c) T-^ (0.25, 0.25, 0)-K. The incident photon energy is tuned to excite the 
Is electron to the peak of the 4p-DOS. The broken and solid lines are the results of the HFA and 
RPA, respectively. The dotted lines represent the experimental data. 10 

The spectra arise from the excitation from the valence states with the e g character to 
the conduction states with also the e g character. Since the unoccupied states are almost 
concentrated in the minority spin states (see Fig. [6]), the electron-hole pair creation occurs 
in the minority spin channel. They form a continuum spectrum. Several peaks are already 
present within the HFA. Some peaks are suppressed and some are enhanced by the RPA. At 
the r point, peaks are found around 4-6 eV and around 8 eV, in good agreement with the 
experiment. In addition, an extra peak is found around the upper edge of the continuum, 
~ 10 eV, which is not confirmed by the experiment (the observation is limited below 10 
eV). With momenta deviating from the T point, intensities around 5-6 eV is particularly 
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enhanced in the entire region of 4-6 eV. Thereby the spectral shape looks like the peak 
position is moving with changing momenta. On the other hand, the peak around 8 eV does 
not change its position with changing momenta. These characteristics are the same along the 
three symmetry lines. The spectral shapes show excellent agreement with the experiment.— 
Note that the sharp peak around 10 eV is strongly suppressed with deviating momenta from 
the T point. 

Below the lower edge of the continuum spectra, we obtain a bound state at 2 eV, as 
shown in Fig. [TUJ This may be compared with the experimental peak at 1.7 eV, but the 
calculated intensity is about two order of magnitude smaller. At present, we do not know 
the real origin for this discrepancy. 



Figure 10: RIXS intensity below the lower edge of the continuum at T point in comparison with the 
experiment.— The calculated intensity is enlarged by multiplying 200 to make the spectra visible. 

V. CONCLUDING REMARKS 

We have analyzed the RIXS spectra in NiO by developing the formalism of Nomura and 
Igarashi. This is based on the Keldysh Green's function formalism and relates the RIXS 
spectra to the density-density correlation function under the Born approximation to the 
core-hole potential. The use of the Born approximation had been justified for the RIXS 
in La2Cu04 by evaluating the multiple scattering effects.— Since NiO has a larger energy 
gap than L^CuO/i, we expect that the Born approximation is better justified in NiO. We 
have introduced the tight-binding model including all the Ni 3d and O 2p orbitals as well 
as the full Coulomb interaction between 3d orbitals, and have calculated the single-particle 
excitation within the HFA. The RIXS spectra are generated by the band-to-band transition 
to screen the Is core hole. We have calculated the RIXS spectra using the single-particle 
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band thus evaluated together with the 4p DOS by the LDA. The electron correlation is 
treated within the RPA. We have obtained several peaks in the continuum spectra. Two 
peaks are found prominent around 5 eV and 8 eV, and the one around 5 eV shifts to the 
higher energy position while another around 8 eV does not move with momenta changing 
to the zone boundary. The calculated spectral shape show quantitative agreement with the 
experiment. The excellent agreement with the experiment suggests that the present scheme 
is useful to analyze the RIXS spectra. As regards the " d-d" excitation, we have obtained 
a bound state below the lower edge of the continuum spectrum. However, the intensity is 
about two order of magnitude smaller than the experimental one. The presence of crystal 
distortion might enhance the intensity of " d-d" excitation. At any rate, the real reason for 
this discrepancy is not known. 
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